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The thermodynamics of the lattice model of intercalation 
of ions in crystals is considered in the mean field approxima- 
tion. Pseudospin formalism is used for the description of 
interaction of electrons with ions and the possibility of hop- 
ping of intercalated ions between different positions is taken 
into account. Phase diagrams are built. It is shown that 
the effective interaction between intercalated ions can lead 
to phase separation or to appearance of modulated phase (it 
depends on filling of the electron energy band). At high val- 
ues of the parameter of ion transfer the ionic subsystem can 
pass to the superfluid-like state. 

1 Introduction 

Theoretical investigation of intercalation of ions in crys- 
tals is an actual problem of modern physics. Metal ox- 
ides as hosts for ion (for example, lithium ions) inser- 
tion are very promising electrode materials. It should be 
noted that theoretical descriptions of such processes in 
most cases were restricted to the numerical ab-initio and 
density-functional calculations. For example, in [1-3] 
quantum-chemical Hartree-Fock and density-functional 
calculations were performed to investigate lithium inter- 
calation in TiC>2 crystal. It was shown that Li is almost 
fully ionized once intercalated (Li looses its valence elec- 
tron) and reconstruction of electron spectrum at inter- 
calation takes place. Thus, ion-electron interaction can 
play a significant role. Another interesting feature of 
such crystals is a shift of the chemical potential at in- 
tercalation into the conduction band. As a result, these 
crystals have metallic conductivity ( [4], for a review, 
see also [5]); before intercalation, such crystals are semi- 



conductors with wide gap. At intercalation of lithium in 
TiC>2, phase separation into Li-poor (Li^o.oiTiC^) and 
Li-rich (Li^o.s-o.eTiCb) phases occurs. This two-phase 
behaviour leads to a constant value of electrochemical 
potential [6,7] (this fact is used when constructing bat- 
teries) . In [8] the Monte Carlo simulation was performed 
to investigate the intercalation using the Hamiltonian 
which included the interaction between ions only. 

In our previous works [9, 10] we have formulated 
the pseudospin-electron model of intercalation and have 
taken into account the ion-electron interaction. It has 
been revealed that the effective attractive interaction 
between ions was formed and the condition of appear- 
ance of phase separation has been established. The 
ion-electron interaction was also considered in [11] at 
the investigation of thermodynamics of S = 1 model 
of intercalation (the model was similar to the known 
Blume-Emery-Griffiths model), but the electron as well 
as ion transfer was not taken into account. It should 
be noted that models of pseudospin-electron model type 
are widely used in physics of the strongly-correlated elec- 
tron systems in recent years. Application of this model 
to high-temperature superconductors allows one to de- 
scribe thermodynamics of anharmonic oxygen ion sub- 
system and explain the appearance of inhomogeneous 
states and the bistability phenomena ( [12]). Such a 
model can also be applied to the description of hydrogen- 
bonded systems. 

In the present paper we deal with a more compli- 
cated model and take into account the possibility of 
the transfer of intercalated ions. The considered model 
corresponds to the hard-core boson approach. Hard- 
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core bosons obey the Pauli spin-1/2 commutation rules. 
Since the original work of Mahan [13], such models were 
applied for description of ionic conductors and calcula- 
tion of their conductivity. Recently one-particle spec- 
trum was investigated in one-dimensional limit [14]. A 
system of hard-core bosons is a particular case of the 
well-known Bose-Hubbard model, which has been inten- 
sively investigated in the last 15 years (see, for exam- 
ple, [15]). The model is of great interest due to the 
experimental realization of optical lattices (see, for in- 
stance, [16]). This model can be directly applied to in- 
vestigate such objects. The Hamiltonian of the Bose- 
Hubbard model includes two terms, one is connected 
with the on-site interaction U between particles and an- 
other term is connected with the particle hopping be- 
tween sites (particles in this model obey the bosonic 
commutation rules). In the limit U — ► oo this model 
reduces to the hard-core boson model. Different theo- 
retical methods were used to study this model: mean- 
field theory [17], random phase approximation [18,19], 
strong coupling approach [20], quantum Monte-Carlo 
method [21]. Recently a bosonic version of dynamical 
mean field theory was formulated [22]. The existence of 
superfluid and Mott-insulator phases is a characteristic 
feature of this model. 



In addition to our previous investigations [9], the aim 
of this work is the study of the ion transfer influence on 
equilibrium states of intercalated ion subsystem. As it 
was shown in [9,10], the effective interaction between 
ions can change (depending on electron band filling) its 
character from repulsion to attraction, leading to the 
charge-ordered modulated phases or phase separation 
into uniform phases with different particle concentra- 
tions, respectively. Ion hopping between local positions 
is unfavorable for the realisation of such phases or phase 
transitions. Besides, ion hopping leads to the appear- 
ance of superfluid type phase. 

We investigate phase transitions in the intercalated 
ion subsystem within the framework of the lattice model 
with ion transfer in the regime of the fixed chemical po- 
tential of the ions and electrons. Electron subsystem is 
described by partially filled one energy band. 



2 The model 

We consider the following model Hamiltonian: 

H = J2n ij SfST+^t ij cfc ja +J2(9S*n ia - f j,n ia )- 

ij ija i<r 

-X>s?. (1) 

i 

Here we introduce the pseudospin variable Sf which 
takes two values; Sf — 1/2 when there is an interca- 
lated ion in a site i and Sf = —1/2 when there is no ion, 
cf a and Ci a are electron creation and annihilation oper- 
ators, respectively. We take into account the possibility 
of ion and electron jumps between sites (the first and 
the second term in (1)) and interaction of electrons with 
ions (g term). The last one is connected with the elec- 
tron band shift at intercalation (such an effect is known, 
for example, for the system Li^TiC^ [3]); [i and h play 
the role of the chemical potentials of electrons and ions. 

It should be noted that we do not consider here the 
direct interaction between ions. In our previous pa- 
per [9], it was shown that ion-electron interaction leads 
to the formation of the effective interaction beweeen ions 
and even at repulsive direct ion-ion interaction the effec- 
tive ion-ion interaction of attractive type can be formed. 
This can lead to the phase transition of the first order 
between uniform phases with jumps of ion and electron 
concentrations. 

The thermodynamics of the model is investigated in 
the mean field approximation (MFA) 

gmS? - g(ni)S? + gn^Sf) - g(m)(S?) 

ns+s- -+ n(s+)s- + nsf(S-) si(s+)(sr), (2) 

here the average ion concentration w — (S z ) + 1/2 is 
introduced; in our approximation (S + ) = (S~) = (S x ), 
(S y ) = 0. The average value (S x ) plays the role of order 
parameter for the case of the superfluid phase (this is 
a phase with condensate of bose-type), and determines 
the concentration of condensed particles. 

Application of the MFA to the strongly correlated 
systems in the limit of a weak one-site correlation makes 
it possible to satisfactorily describe their properties, 
when there is no correlational splitting of the electron 
band. This approximation is an analogy to the virtual 
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crystal approximation, which is often used for mixed 
systems. To go beyond the MFA one can use more 
complicated approximations, for instance, the coherent 
potential-like approximations. Besides, in the case of 
the Bose-Hubbard model the kinetic energy term is of- 
ten considered within the mean field approach. This ap- 
proximation is well known to give a reasonable estimate 
of the critical on-site repulsion at which Mott-insulator 
- superfluid phase transition occurs [15,17]. 

In [23, 24] the case O = was considered. It was 
shown that if the chemical potential is near the band 
center, the double modulation phase is realised in the 
system, while in the case when the chemical potential 
is close to the band edges, the phase transition be- 
tween uniform phases occurs. At intermediate values of 
the chemical potential the incommensurate modulated 
phase appears. In the present investigation we restrict 
ourselves to the cases of double modulation phase and 
uniform phase. 

The Hamiltonian of the model in the MFA is as fol- 
lows: 



H MFA = ^(g^ - M ) ni , 



]T(. 9 n Q - h)S z ia + 



- E *t c t*, c iP<> + E 2^ /3 <^>^ - g E n aVa - 

ia,j{3 otfdi ia 

-Nn(Sf)(s$), 



(3) 



here we consider two sublattices: (J2 a n iaa) = 
n ai {Sfa) = Va', ot — 1,2 is a sublattice index, i 
is an unit cell index, N is the number of lattice sites, 

n = n 12 = n 21 = v. n} 2 ; n aa = t aa = o. 

This Hamiltonian can be diagonalized. We pass to k- 
representation and perform the unitary transformation 
in the pseudospin subspace 



H 



" ' A = E(^ kQ _ M)"kacr - E ^a°fc 



acrk 



g^L(n im +n 2m )-Nn{Sl){S%) 



(4) 



A kQ = g — ^ — + (-1) \ (g 



Cki<T = Ckio- cos 6 + Ck2<T sin < 



Ck2a = -Ckia sin 6 + C k2(T COS ( 

tk 



sin 26 



S? a = °ia cos 6 a + af a sin 9 a , 



s f a = °fa cos 9 a - o z ia sin 6 a , sin 6 a = 



2fi(S|> 



\ a = ^(gn a -h) 2 + (2f!(^)) 2 , a?0. 

The doubling of unit cell leads to the splitting in the 
electron spectrum. Two subbands are separated by the 
gap g\rji — 772 ] - The electron band changes its position 
at intercalation. 

Using MFA, we can calculate the mean values of both 
the electron and ion concentrations: 

1 ^ 1 + cos 20 ^q-c 



N ^ 2 

ker 



(e^" +1)- 1 - 



1 - cos 26 , , 



(5) 



Va = — ^ tanh(— ),{S%) = -^tanh(— ). 



To find the thermodynamically stable states we also 
have to calculate the grand canonical potential 

■W = -m E ln (( e ^" + !)" )( e ^" + !)" ))- 

2 kcr 



-Tln(4cosh(^- cosh(^2.)) - g{n 1 r]i + n 2 77 2 )- 



2fi<Sf)<Sf>. 



(6) 



The absolute minima of the <&-function determine the 
equilibrium states. 
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3 Results 

The semielliptic density of states, p(e) = j^VW 2 — e 2 , 
—W<e<W, where W is a half width of the elec- 
tron band, was used (W is chosen as energy unit; in 
our calculations we put W = 1). Using this density of 
states, we perform summation over k in the equations of 
self-consistency J5]) and in the expression for the grand 
canonical potential (jfij. 

As noted above, the stable states are obtained using 
the condition of minimum of the function $. In Fig.l 
and Fig. 2 the (h — f2) phase diagrams are shown for the 
cases [i = (at the center of the band) and /i = — 0.7W 
(near the lower band edge). 



two uniform phases with jumps of the average ion and 
electron concentration (and the phase separation in the 
regime of the fixed concentration, see [9, 10, 23, 24] for 
more details) takes plase (dashed line in Fig. 2). 

It is easy to see that at high values of the param- 
eter of ion transfer Q the only possible phases are the 
uniform phases with (S x ) ^ and (S x ) = 0. The phase 
with (S x ) 7^ appears due to the presence of ion hopping 
between sites; this phase is an analogy to a superfmid 
phase in the systems of hard-core bosons and can cor- 
respond to the state with high mobility of intercalated 
ions. In the case ft < 0, one finds (Sf) = (5f), while in 
the case > 0, one finds (Sf) = -(Sf). 
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Fig.l Phase diagram in the (h — fi) plane at /i = 0. The parameter 
values are: g = -0.4, W = 1,T = 0.03. (1),(3) - uniform phase 
with (S*> + ((1): (Sf } = (Sf >, (3): (Sf ) = -(Sf »; (2)-uniform 
phase with (S x ) = 0; (4)-modulated phase ((a): (Sf 2 ) = 0, (b): 
(Sf 2 ) 7^ 0). The solid line denotes the phase transitions of the 
second order and the dotted line denotes the phase transition of 
the first order. 

In the case fj, = and at small values of ion hop- 
ping parameter 12 the system undergoes the phase tran- 
sition of the first order from the uniform to modulated 
phase (in the modulated phase n% ^n,2,T)i ^ 772) at the 
change of the chemical potential of the ions (dotted line 
in Fig.l). In the case fi = — 0.7W and at small val- 
ues of the phase transition of the first order between 
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Fig. 2 Phase diagram in the (h — O) plane at /1 = —0.7. The pa- 
rameter values are: g = —0.4, W = 1, T = 0.03. The notations are 
the same as in Fig.l. The solid line denotes the phase transition of 
the second order and the dashed line denotes the phase transition 
of the first order. 

In Fig. 3 and Fig. 4 the [h — fx) phase diagrams are 
shown for the cases O = and £1 = 0.25. Dotted line 
denotes the first order phase transition between uniform 
and modulated phases, dashed line denotes the first or- 
der phase transition between uniform phases. 
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Fig. 3 Phase diagram in the (h — fi) plane at Q = 0. The parame- 
ter values are: g = —0.4, W = 1, T = 0.03. The notations are the 
same as in Figs. 1,2. 



This curve ends in critical point at some value of tem- 
perature T cr . At high values of ft this phase transition 
disappears (see Figs. 2, 6). The existence of such phase 
transitions (in the regime of the fixed concentration it 
corresponds to the phase separation into phases with 
different concentration of ions) is in accordance with ex- 
perimental data for intercalated crystals, where the ap- 
pearance of poor and rich ion concentration phases was 
observed (see, for example, [6]). The presence of mod- 
ulated phase in intercalated crystals is also indicated in 
experiments (for review, see [5]). 

The phase transition from the uniform to modulated 
phase can be of the second or the first order, this is il- 
lustrated in Figs. 7,8 for the cases O = and ft = 0.3. 
We should draw attention to the fact, that with increas- 
ing temperature the first order transition will transform 
into the second order one and then will disappear (as it 
is shown in Figs. 7, 8). 




The next four diagrams, which are shown in Figs. 
5-8 are the (T— h) phase diagrams. In Fig. 5 the first or- 
der phase transition curve between two uniform phases 
with jumps of ion and electron concentrations is shown. 



Fig. 5 Phase diagram in the (T — h) plane at U = 0. The parame- 
ter values are: g = —0.4, W = 1, /X = —0.7. The notations are the 
same as in Figs. 1,2. 
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Fig. 6 Phase diagram in the (T — h) plane at Q = 0.2. The param- 
eter values are: g = — 0.4, W = 1,/i = —0.7. The notations are 
the same as in Figs. 1,2. 
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Fig. 8 Phase diagram in the (T — h) plane at f! = 0.3. The param- 
eter values are: g = — 0.4, W = 1, /i = 0. The notations are the 
same as in Figs. 1,2. 

4 Conclusions 
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Fig. 7 Phase diagram in the (T — h) plane at Q = 0. The parameter 
values are: g = — 0.4, W = = 0. The notations are the same 
as in Figs. 1,2. 



In this work, the pseudospin-electron model of ion inter- 
calation in crystals has been formulated. The model can 
be applied for the description of the thermodynamics of 
such a process in materials with electron bands where 
the band filling has metallic or semimetallic character 
(in particular, compounds of transition metals such as 
Ti02 or other similar systems with narrow conduction 
bands). The thermodynamics of the model has been 
investigated in the mean-field approximation. The effec- 
tive interaction between intercalated ions is formed due 
to their interaction with electron subsystem. Such an 
interaction is attractive or repulsive depending on the 
filling of electron band (in the first case the chemical 
potential of electrons should be close to the band edge; 
the second case is realised near half filling) . The appear- 
ance of modulated phase or phase transitions of the first 
order with jumps of ion and electron concentrations (in 
the regime of the fixed concentrations it corresponds to 
the phase separation) has been established. 

Increase of the ion transfer parameter leads to the 
disappearance of both of modulated phase and phase 
transition with jumps of the ion and electron concentra- 
tions. Besides, the new phase with (S x ) ^ appears due 
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to ion hopping between sites; this phase is an analogy 
to superfluid phase in the systems of hard-core bosons 
or superionic phase in the crystalline ionic conductors 
(phase with high mobility of ions). Such a phase can 
exist at intermediate values of the chemical potential of 
intercalated ions and the transition to this phase is of 
the second order. To investigate this phase in detail we 
should examine the behaviour of the conductivity and 
other characteristics of the system. This is the task for 
future investigations. 

References 

[1] A. Stashans, S. Lunell, R. Bergstrom, A. Hagfeldt, 
S.-E. Lindquist, Phys. Rev. B 53, 159 (1996). 

[2] M.V. Koudriachova, N.M. Harrison., W. de Leeuw, 
Phys. Rev. Lett. 86, 1275 (2001). 

[3] M.V. Koudriachova, N.M. Harrison, W. de Leeuw, 
Phys. Rev. B 65, 235423 (2002). 

[4] V. Luca, B. Hunter, B. Moubaraki, K.S. Murray, 
Chem. Mater. 13, 796 (2001). 

[5] I.I.Grygorchak, Phys. Chem. Sol. State 2, 7 (2001) 
(in Ukrainian). 

[6] M. Wagemaker, R. van de Krol, A. P.M. Kentgens, 
A. A. van Well, F.M. Mulder, J. Am. Chem. Soc. 
123, 11454 (2001). 

[7] M. Wagemaker, G.J. Kearley, A. A. van Well, H. 
Mutka, F.M. Mulder, J. Am. Chem. Soc. 125, 840 
(2003). 

[8] K.-N. Jung, S.-I. Pyun, S.W. Kim, J. Power 
Sources 119-121, 637 (2003). 

[9] T.S. Mysakovych, I.V. Stasyuk, J. Phys. Studies 
11, 195 (2007). 

[10] T.S. Mysakovych, V.O. Krasnov, I.V. Stasyuk, 
Condens. Matter Phys. 11, 663 (2008). 



[11] I.V. Stasyuk, Yu.I. Dublenych, Phys. Rev. B 72, 
224209 (2005). 

[12] I.V. Stasyuk, Order, disorder and criticality . Ad- 
vanced problems of phase transition theory. Vol- 
ume 2, edited by Yu. Holovatch (World Scientific, 
Singapore, 2007), p. 231. 

[13] G.D. Mahan, Phys. Rev. B 14, 780 (1976). 

[14] I.V. Stasyuk, I.R. Dulepa, Condens. Matter Phys. 
10, 259 (2007). 

[15] N. Dupuis, K. Sengupta, Physica B: Physics of 
Condensed Matter 404, 517 (2008). 

[16] M. Greiner, O. Mandel, T. Esslinger, T.W. 
Hansen, I. Bloch, Nature 415, 39 (2002). 

[17] K. Sheshadri, H.R. Krishnamurthy, R. Pandit, 
TV. Ramakrishnan, Europhys. Lett. 22, 257 

(1993) . 

[18] S. Konabe, T. Nikuni, M. Nakamura, Phys. Rev. 
A 73 , 033621 (2006). 

[19] Y. Ohashi, M. Kitauri, H. Matsumoto, Phys. Rev. 
A 73, 033617 (2006). 

[20] J.K. Freericks, H. Monien, Europhys. Lett. 26, 545 

(1994) . 

[21] B. Capogrosso-Sansone, S. G. Soyler, N. 
Prokof'ev, B. Svistunov, Phys. Rev. A 77, 
015602 (2008). 

[22] K. Byczuk, D. Vollhardt, Phys. Rev. B 77, 235106 
(2008). 

[23] I.V. Stasyuk, T.S. Mysakovych, J. Phys. Studies 
5, 268 (2001). 

[24] I.V. Stasyuk, T.S. Mysakovych, Condens. Matter 
Phys. 5, 473 (2002). 



7 



